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With the help of variational perturbation theory we con- 
tinue the renormalization constants 4 -theories in 4 — e dimen- 
sions to strong bare couplings go and find their power behav- 
ior in go, thereby determining all critical exponents without 
renormalization group techniques. 

1. In a recent paper [El we have shown that there exists 
a simple way of extracting the strong-coupling properties 
of a </3 4 -theory from perturbation expansions. In partic- 
ular, we were able to find power behavior of the renor- 
malization constants in the limit of large couplings, and 
from this all critical exponents, of the system. By us- 
ing the known expansion coefficients of the renormaliza- 
tion constants in three dimensions up to six loops we de- 
rived extremely accurate values for the critical exponents. 
The method is a systematic extension of the Feynman- 
Kleinert variational approximation to path integrals j|] 
to arbitrary orders For an anharmonic oscillator, 
the derived variational perturbation expansions converge 
uniformly and exponentially fast, like e - constxN /3 m the 
order N of the approximation ^,|). The same type of 
convergence seems to exist also for the </> 4 -theory if the 
power 1/3 is replaced by 1 — to, where u> is the criti- 
cal exponent governing the approach to scaling [Q . This 
exponent plays a crucial role in the development of the 
theory. 

2. Variational perturbation expansions have the im- 
portant property of possessing a good strong-coupling 
limit as was first shown for the harmonic oscillator |^|,^| . 
The speed of convergence turned out to be governed 
by the convergence radius of the strong-coupling expan- 
sion dJ^]. The good strong-coupling properties have en- 
abled us to set up a simple algorithm for deriving uni- 
formly convergent approximations to functions of which 
one knows a few Taylor coefficients and an important 
scaling property: they approach a constant value with 
a given inverse power of the variable. The renormalized 
coupling constant g of a </5 4 -theory has precisely this prop- 
erty as a function of the bare coupling constant go- In 
D = 4 — e dimensions, it approaches a constant value g* 
for increasing bare coupling constant go like 
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where g* is the infrared-stable fixed point and lu is called 
the critical exponent of the approach to scaling. 

The purpose of this paper is to point out that the the- 
ory developed in jjj] for a three-dimensional 4 -theory can 
easily be applied in D = 4 — e dimensions with beautiful 
results at the two-loop level. 

3. Let us briefly recall the relevant formulas. Consider 
a function f(go) for which we know TV expansion terms, 

/jvCffo) — Y,n=o a n9o , an d the fact that it approaches a 
constant value /* in the form of an inverse power series 

/m(5o) = J2m=o b -m(9o 2/q ) m witn a nnite convergence 
radius g s (simple examples were treated in [fiof ). Then 
the TVth approximation to the value /* is obtained from 
the formula 
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where the expression in brackets has to be optimized 
in the variational parameter go- The optimum is the 
smoothest among all real extrema. If there are no real 
extrema, the turning points serve the same purpose. 

The derivation of this expression is simple: We replace 
.go in /zv(flo) trivially by g = g /n q with k = 1. Then we 
rewrite, again trivially, n~ q as (K 2 + k 2 — K 2 )~ q / 2 with 
an arbitrary parameter K. Each term is now expanded in 
powers of r = (k 2 — K 2 ) / K 2 assuming r to be of the order 
go . Taking the limit go — > 00 at a fixed ratio go = go/ K q , 
so that K — > 00 like g\J q and r — > —1, we obtain (||). 
Since the final result to all orders cannot depend on the 
arbitrary parameter K, we expect the best result to any 
finite order to be optimal at an extremal value of K, i.e., 
of g . 

The strong-coupling approach to the limiting value 
r = — 1 + k 2 /K 2 = — 1 + 0(gQ 2 ^ q ) implies the leading 
correction to /jy to be of the order of g Q 2 ^ q . Applica- 
tion of the theory to a function with the strong-coupling 
behavior (l|) requires therefore a parameter q — ltjui in 
formula (p). 

For N — 2 and 3 one can give analytic expressions for 
the strong-coupling limits (|^). Setting p = 1 + q/2 = 
1 + e/ui, we find for N = 2 
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For vV = 3, we obtain from the extrema 
/s = °Pt K + i a iP(p + l)<7o + a 2 {2p - 1)^ + a 3 gl] 
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where r = yl — 30103/02 and ai = \a\p{p + 1) and 
d 2 = a 2 (2p — 1). The positive square root must be taken 
to connect 3J smoothly to g% in the limit of a vanishing 
coefficient of 35*. If the square root is imaginary, the 
optimum is given by the unique turning point, leading 
once more to (||) but with r = 0. 

The parameter p = 1 + e/u; can be determined from 
the expansion coefficients of a function F(go) as follows. 
Assuming F(go) to be constant F* in the strong-coupling 
limit, the logarithmic derivative /(go) = 9oF / (go)/F(go) 
must vanish at g = 00. If F(g ) starts out as A +Aig + 
. . . or Aigo + A 2 g^ + . . . , the logarithmic derivative is 



/(So) = A[ go + (2A' 2 - A[ 2 )g 2 

+ K 3 - 3A[A' 2 + 3A 3 ) 5o 3 

where A' i = Ai/A , or 

f(g ) = 1 + A 2 g + (2i 3 - A\)gl 
+ (i 3 - 3A 2 A 3 + 3A 4 )g 3 - 
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where Ai — Ai/A\. The expansion coefficients on the 
right-hand sides are then inserted into (|J) or (Q) , and the 
left-hand sides have to vanish to ensure that F(go) — > F*. 

If the approach F(g ) — * F* is of the type (0), the 
function 
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must have the strong-coupling limit 
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4. These formulas are now applied to the renormaliza- 
tion constants of the c6 4 -theory in D = 4 — e dimensions 
with the bare euclidean action 



.4: : ld D X {\[dU^f+\<^) + (^) 2 ^ [0O(*)] : 



The field c6o(x) is an n-dimensional vector, and the action 
is 0(n)-symmetric in this vector space. The Ising model 
corresponds ton = 1 , the critical behavior of percolation 
is described by n — 0, superfluid phase transitions by 
n = 2, and classical Heisenberg magnetic systems by n = 
3. 

By calculating the Feynman integrals regularized via 
an expansion in e = 4 — D with the help of an arbitrary 
mass scale p, one obtains renormalized values of mass, 
coupling constant, and field related to the bare input 
quantities by renormalization constants Z$, Z m , Z g : 



m\ = m 2 Z m Z, x , \ Q = \Z g Z, 2 , c6 = zV 2 . (10) 



Up to two loops, perturbation theory yields the follow- 
ing expansions in powers of the dimensionless reduced 
coupling constant g = X/p e : 
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We now set the scale parameter p equal to m and con- 
sider all quantities as functions of go — A/m e . In order 
to describe second-order phase transitions, we let m 2 , go 
to zero like r =const x (T — T c ) as the temperature T ap- 
proaches the critical temperature T c . Then also m 2 will 
go to zero, and thus go to infinity. Assuming the theory 
to scale as suggested by experiments, we now determine 
the value of the renormalized coupling constant g in the 
strong-coupling limit go — > 00, and the approach to it, as- 
suming the behavior (|l|). For this we apply formula (p) 



to 3(30) to find g* , and use the vanishing of (^]), or (|) 
with (|7p at strong couplings to determine w. Under the 
scaling assumption, the ratios rn 2 /rn^ and c6 2 /c6q have 
the limiting power behavior for small m: 
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The powers can therefore be calculated from the strong- 
coupling limits of the logarithmic derivatives 
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From @ and Q find the expansions 
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When approaching the second-order phase transitions, 
(9) where rrip vanishes like r = (T — T c ), m 2 vanishes with a 
different power of r. This power is obtained from the first 
equation in (|IJ) which shows that m oc r 1 / ( 2- f"») . Exper- 
iments observe that the coherence length of fluctuations 
£ = 1/m increases near T c like t"^, so that we derive for 
the critical exponent v a value 1/(2 — rj m ). Similarly we 
see from the first equation in ([lj) that the scaling dimen- 
sion D/2 — 1 of the free field 4>o for T — > T c is changed, 
in the strong-coupling limit 30 — * o°i to _D/2 — 1 + r//2, 
the number 77 being the so-called anomalous dimension of 
the field. This implies a change in the large-distance be- 
havior of the correlation functions ((f)(x)<f>(0)) at T c from 
the free-field behavior r~ D+2 to r - £) + 2 -'). The magnetic 



2 



susceptibility is determined by the integrated correlation 
function ((f>o(x)<j>o(0)}. At zero coupling constant go, this 
is proportional to 1/m 2 , oc t -1 , which is changed by fluc- 
tuations to m -2 (f>Q I '(j) 2 . This has a temperature behavior 
TO -(2-'f) = r -" y (2-»)) = t~ 7 , which defines the critical ex- 
ponent 7 = v{2 — rj) observable in magnetic experiments. 
Using v — 1/(2 — r) m ) and the expansions ([i"6|), (p7|), we 
obtain for 7(50) the perturbation expansion up to second 
order in go : 
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All calculations in this note will be restricted to the 
two loop level, which will be sufficient to demonstrate 
the power and beauty of the new strong-coupling theory 
with analytical results. 

5. We begin by calculating the critical exponent u> 
from the requirement that g(go) has a constant strong- 
coupling limit, implying the vanishing of (|^) for go — > 
00. From the expansion ( pi] ) we obtain a logarithmic 
derivative (^) up to the term g^, so that Eq. ([|) can be 
used to find the scaling condition 
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This gives 
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Since u> must be greater than zero, only the positive 
square root is physical. With the explicit coefficients 
Ai,A2,As of expansion (11), this becomes 



3n + 14 
p = 2W1 + 3 7 —re. 
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The associated critical exponent to = e/(p — 1) is plotted 
in Fig. 1. It has the e-expansion 
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which is also shown in Fig. 1, and agrees with the first 
two terms obtained from renormalization group calcula- 
tions @. 

From Eqs. (||), (f7|), and (||) we obtain for the critical 
exponent w a further equation 




with the positive sign of the square root ensuring a pos- 
itive to. Inserting the coefficients of (fl]), this becomes 



(25) 



The associated critical exponent u> = e/(p — 1) has the 
same e-expansion ( |22| ) as the previous approximation 
(pl|). The full approximations based on ( f25| ) is indis- 
tinguishable from the earlier one in the plot of Fig. 1. 

Having determined u>, we can now calculate g* . Insert- 
ing the first two coefficients of the expansion (jl^) into (|^) 
we obtain 
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Inserting (pi|), this yields 
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which is precisely the well-known e-expansion of g* in 
renormalization group calculations up to the second or- 
der. Including the next coefficient, we can use for- 
mula (^) to calculate the next approximation g 3 . At 
e = 1, the square root turns out to be imaginary, so 
that it has to be omitted (corresponding to the turn- 
ing point as optimum). The resulting curve lies slightly 
(« 8%) above the curve (p7|), i.e., represents a worse 
approximation than (|2~7j). Indeed, the e 3 -term in g% 
is 81(3n + 14) 2 /8(n + 8) 5 and disagrees in sign with 
the exact term e 3 [3 (-33n 3 + 110n 2 + 1760n + 4544) /8 
-36C(3)(n + 8) (5n + 22)]/(n + 8) 5 , which we would find 
by calculating p from an expansion ([ll]) with one more 
power in go- 

We now turn to the critical exponent v. Taking the 
expansion ( |l6| ) to go — > 00, we obtain from formula (j^) 
the limiting value 
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The corresponding v — 1/(2 — rj m ) is plotted in Fig. 1. 
With the approximation (|2lj) for p we find for v the e- 
expansion 
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which is solved by 
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which is also shown in Fir. 1, and agrees with renormal- 
ization group results to this order. 

As a third independent critical exponent we calculate 
7 = (2 — n)/(2 — rjm) by inserting the coefficients of the 
expansion (Jig) into formula (H), which yields 
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plotted in Fig. 1. This has an e-expansion 
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shown again in Fig. 1, and agreeing with renormalization 
group results to this order. The full approximation is 
plotted in Fig. 1. The critical exponent r\ — 2 — ~f/v has 
the e-expansion rj = (n + 2)e 2 /2(n + 8) 2 + . . . . 

6. We conclude that variational strong-coupling the- 
ory can easily be applied to 4 -theories in D = 4 — e 
dimensions and yields resummed expressions for the e- 
dependence of all critical exponents. Their e-expansions 
agree with those obtained form renormalization group 
calculations. The calculations here are supposed to il- 
lustrate the new calculational procedure and do not yet 
give good results for the critical exponent uj. In order to 
achieve high accuracy, we shall have to go to five loops 
and incorporate knowledge of the large-order behavior of 
the expansion coefficients. 

Recent results by A. Pelissetto and E. Vicari (Univer- 
sity of Pisa preprint IFUP-TH 52/97) on the renormal- 
ization in 0(n)-models via 1/n expansions will lead to 
further improvements of the theory 
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FIG. 1. For the Ising universality class (n = 1), the first 
figure shows the renormalized coupling at infinite bare cou- 
pling as a function of e = 4 — D calculated via variational 
perturbation theory from the first two perturbative expan- 
sion terms. The curve coincides with the e-expansion up 
to order e 2 . The dashed curve indicates the linear term. 
The other figures show the critical exponents uj, v, and 7. 
Dashed curves indicate linear and quadratic e-expansions. 
The dots mark presently accepted values of g* ~ 0.48 ±0.003, 
uj » 0.802 ± 0.003, v = 0.630 ± 0.002, and 7 = 1.241 ± 0.004 
obtained from six-loop calculations M. 



5 



